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The aim of the paper is to understand how the inclusion of more and more time-scales into a 
stochastic stationary Markovian process affects its conditional probability. To this end, we consider 
two Gaussian processes: (i) a short-range correlated process with an infinite set of time-scales 
bounded from below, and {ii) a power-law correlated process with an infinite and unbounded set of 
time-scales. For these processes we investigate the equal position conditional probability P{x, t\x, 0) 
and the mean First Passage Time Tx{A.). The function P{x,t\x,Q) can be considered as a proxy 
of the persistence, i.e. the fact that when a process reaches a position x then it spends some time 
around that position value. The mean First Passage Time can be considered as a proxy of how fast 
is the process in reaching a position at distance A starting from position x. 

In the first investigation we show that the more time-scales the process includes, the larger the 
persistence. Specifically, we show that the power-law correlated process shows a slow power-law 
decay of P{x,t\x,Q) to the stationary pdf. By contrast, the short range correlated process shows a 
decay dominated by an exponential cut-off. Moreover, we also show that the existence of an infinite 
and unbonded set of time-scales is a necessary and not sufficient condition for observing a slow 
power-law decay of P{x, t\x, 0). In fact, in the context of stationary Markovian processes such form 
of persistence seems to be associated with the existence of an algebraic decay of the autocorrelation 
function. In the second investigation, we show that for large values of A the more time-scales the 
process includes, the larger the mean First Passage Time, i.e. the slowest the process. On the other 
hand, for small values of A, the more time-scales the process includes, the smaller the mean First 
Passage Time, i.e. when a process statistically spends more time in a given position the likelihood 
that it reached nearby positions by chance is also enhanced. 

PACS numbers: 02.50.Ey, 05.10. Gg, 05.40.-a, 02.50.Ga, 89. 75. Da 



I. INTRODUCTION 

Stochastic processes are used to model a great variety 
of systems in disciplines as disparate as physics [iHS], 
genomics [3 [5], finance [nUTU], cHmatology [TT] and social 
sciences [T^ . 

Among others, Markov processes play a central role 
in the modeling of natural phenomena due to their sim- 
plicity. In fact, a Markov process is fully determined 
by the knowledge of its probability density function 
(pdf) W{x,t) and its conditional transition probability 
P{xn+i, tn+i\xn, tn)- When a stationary Markovian pro- 
cess is continuous both in space and time, the time evolu- 
tion of these distributions is described by a Fokker-Planck 
(FP) equation [HIS]. In this paper we will only consider 
stationary Markovian processes that can be described by 
using a FP equation. 

Stochastic processes can be considered from the point 
of view of their correlation properties. Under this re- 
spect, random variables are divided in short-range and 
long-range correlated variables. Short-range correlated 
variables are characterized by a finite mean of time-scales 
of the process whereas a similar mean time-scale does not 
exist for long-range correlated variables [TBHTT] . 

Several stationary Markovian processes are short-range 
correlated. In fact, the paradigmatic stationary Marko- 
vian process is the Ornstein-Uhlenbeck (OU) one [18], 
whose autocorrelation function is the exponential func- 
tion p(t) = e~'^/'^ where T is the time-scale of the pro- 



cess. However, generally speaking a stationary Marko- 
vian process can be multi-scale, i.e. it may admit ei- 
ther a discrete or a continuum set of time-scales [T21- 
[25] . In general, the autocovariance function R{t) = 
{x{t)x{t + r)) of an additive stationary Markovian pro- 
cess admitting a Fokker-Planck equation can be written 
as R{t) = En=o exp(-AT) + /_^'^c^exp(-Ar), where 
Aq — 0, Ai,...,Ap and ]Ac,+oo[ (Ac ^ Xp) constitute the 
discrete and continuum parts of the eigenvalue spectrum 
of the FP equation. The quantities Cn, c\ are appropri- 
ate weights computed in terms of the eigenfunctions [3J . 
This picture indicates that the more time-scales included 
in the process, the slower the decay of the autocorrela- 
tion function. A generalization to non additive processes 
that are obtained by performing appropriate coordinate 
transformations of the additive ones has been considered 
inRef. [M]- 

Our main interest is to investigate what are the deter- 
minants of a long-range correlated stationary Markovian 
process, at the level of its probability distributions. In 
Ref. [51] we considered stationary Markovian stochas- 
tic processes with a power-law decaying autocorrelation 
function and a stationary pdf with Gaussian or expo- 
nential tails. Those results suggested that the station- 
ary pdf does not bear any special feature with respect 
to the correlation properties of the process. We there- 
fore want to investigate how the inclusion of more and 
more time-scales into stationary Markovian stochastic 
processes affects their conditional probability, which is 
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the other probability distribution that fully characterizes 
the process. 



II. PROCESSES WITH GAUSSIAN 
STATIONARY PDF 



A. A short-range correlated process 



To this end we will consider two stationary Markovian 
processes with Gaussian pdf and autocorrelation func- 
tion not described by a single exponential. In one case 
the process is short-range correlated and in the other 
case the autocorrelation function decays like a power- 
law. We preliminarily show that, as suggested by the 
Doob theorem (25j, this is only possible if one admits a 
conditional probability different from Gaussian for both 
processes. We then investigate the conditional probabil- 
ity by following two complementary approaches. First we 
consider the conditional probability at equal spatial posi- 
tions P{x, t\x, 0) as a proxy of how persistent is a process 
in a given position. We show that the higher the number 
of time-scales incorporated in the stochastic process the 
higher the persistence. 



Let us consider the stochastic process described by the 
following Langcvin equation [3]: 



x{t) = -Kxit)) +DT{t) 
+fc if X < 0, 

h{x) ■ 

-fc if x> 0. 



(1) 



where fc is a real constant and T{t) is a (5-correlated Gaus- 
sian noise term. The diffusion coefficient D will be set 
to unity hereafter. The stationary distribution of this 
process is: 



W,{x) 



■ exp(— fc|a;|) 



(2) 



By using the methodology of eigenfunction expansion [31 
H] it is possible to prove that the autocovariance function 
R{t) — {x{t)x{0)) of the above process is: 



We then consider a quantity related to the conditional 
probability, i.e. the mean First Passage Time T^iA) [3]. 
It can be considered as a proxy of how fast is the pro- 
cess in reaching a position at distance A from the initial 
position x. We show that the higher the number of time- 
scales incorporated in the stochastic process the slower 
the process is. On the other hand, for small values of 
A, the more time-scales the process includes, the smaller 
the mean First Passage Time, i.e. when a process sta- 
tistically spends more time in a given position the likeli- 
hood that it reached nearby positions by chance is also 
enhanced. 



The paper is organized as follows: in section |TT] we will 
show two examples of processes that are Gaussian dis- 
tributed and whose autocorrelation function is not ex- 
ponential. In one case, see section |II A( the process is 
short-range correlated. In the second case, see section 
|IIB[ the process can be long-range correlated with an 
appropriate choice of the relevant parameters. In any 
case the process is always power-law correlated. In sec- 
tion [TiT] we consider the conditional probability at equal 
spatial position P(x,t\x,0) and show that it decays like 
a power-law decays for the algebraically correlated pro- 
cess. We also show that such peculiar behaviour is not 
shared by other processes that include an infinite and un- 
bounded set of time-scales. In section|IV|we consider the 
mean First Passage Time (L) and show that the higher 
the number of time-scales incorporated in the stochastic 
process the slower the process is. Finally, in section |V] 
we draw our conclusions. 



^ ^(1 - 2^ + 4^' + l^') (1 - Erf(V^) + (3) 



^^(2r-l)(3 + 2r)exp(-r) 



fc2 

T = —t 
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which behaves like a power-law with an exponential trun- 
cation for large time lags: Rs{t) ~ exp{—^t)t~^/^ as 
t — !■ oo. 

Let us now consider the coordinate transformation: 



-k\x\l 



(4) 



By using the Ito lemma, one can show that, starting from 
the process of Eq. ([T]) , in the coordinate space y — fs (x) 
one gets a multiplicative stochastic process whose sta- 
tionary pdf is given by: 



Wg{y) 



1 . 1 2^ 

'2tts 2s 



(5) 



with s an additional arbitrary parameter. The corre- 
sponding drift and diffusion coefficients are given by: 



i7(y) = ifc2e^/V(2^) (l- Erf (12/1/72^ 



(6) 



X (^TTp/Z^^s) (^i_Erf(|y|/%/2l))-e2%/2"7rs; 
Giy) = fcy^e^^/(2^) (l-Erf(|y|/V2;)) (7) 

where e = -1-1 when y > and e = —1 when y < 0. For 
large values of y the drift and diffusion coefficient behave 
as: 



H{y) cx -k^s 



1 



G{y) cx +ks 



1 



-oo. 
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FIG. 1: The figure shows the drift coefficient H{y) (top panel) 
and diffusion coefficient G{y) (bottom panel) of the process 
defined by Eq. ([TJ with the coordinate transformation of Eq. 
Q for the case when k = 2.0 and s = 1.0, which corresponds 
to a variance {y{tY) = 0.5. 



In Fig. [Tjwe show the drift coefficient H{y) (top panel) 
and diffusion coefficient G{y) (bottom panel) for the case 
when k = 2.0 and s — 1.0, which corresponds to a vari- 
ance {x{tf) = 0.5 and {y{t)'^) = 1.0. The diffusion co- 
efficient G{y) is continuous in y — Q although its first 
derivative is discontinuous. The drift coefficient Hiii) 
suffers a discontinuity in y = 0. 

According to the eigenfunction expansion methodology 
[31 S] , the autocorrelation function of the process defined 
by Eq. ([T]) with the coordinate transformation of Eq. Q 
is given by p,{t) = {{y{t + T)y{t)) - {yit))^) / {{y\t)) - 
(?/(t))^) where: 



{y{t)y{t + T)) 



(8) 



Cx = 



+ 00 



dx fs{x) ipaix) tpxix) 



where "00(2;) and ip\{x) are the eigenfunctions of the 
Schrodinger equation with potential Vs{x)= k'^/A—k S{x) 
associated to the stochastic process of Eq. ([I]). Eq. Q 
can be used to numerically obtain the autocorrelation of 
the new process defined by Eq. ([ij with the coordinate 
transformation of Eq. Q . 

In the top panels of Fig. [2] we show the results of nu- 
merical simulations of the autocorrelation function per- 
formed for the case when s = 1.0, and k = 2.0 (left), 
k = 0.1 (right). The solid (red) line shows the theo- 
retical predictions obtained from Eq. ([s]), while the open 
circles show the results of the numerical simulations. The 
bottom panels of Fig. [2] show the stationary pdf of the 
process when s = 1.0, and k = 2.0 (left), k = 0.1 (right). 
Again the solid (red) line shows the theoretical predic- 



FIG. 2: (Color online) The figure shows time-average numer- 
ical simulations performed according to Eq. ( 10 1 for the case 
when when s = 1.0, and k = 2.0 (left), k — 0.1 (right). In the 
top panel we show the results for the autocorrelation function. 
The solid (red) line shows the theoretical predictions obtained 
from Eq. Q, while the open circles show the result of the nu- 
merical simulations. The bottom panel shows the stationary 
pdf of the process. Again the solid (red) line shows the theo- 
retical prediction of Eq. ([5|, while the open circles show the 
result of the numerical simulations. The parameters of the 
numerical simulations are given in the text. 



tion of Eq. ^ , while the open circles show the result of 
the numerical simulations. 

Those shown in Fig. [2] are time-average numerical sim- 
ulations performed according to the relation: 



p{t) = — / dt x^.{t)x^.{t + t) 
^ Jo 



(9) 



where T is the length of the simulated time-series and 
X* (t) is one realization of the process. Indeed, in order 
to improve the statistical reliability of our numerical sim- 
ulations, in the region r > 1 we have also averaged over 
a number M of different realizations of the process: 



Pt{t) 



dtXj(t)x,{t + T) 



(10) 



The data shown in the figure are the mean and the stan- 
dard deviations of the M autocorrelation values com- 
puted in each iteration for each time lag. The param- 
eters of the numerical simulations of the autocorrelation 
function are: M = 1 in the region r < 1, Af = 100 in 
the region t e [1, 10] and M = 5 in the region r > 10 
(not shown in the panel) when s = 1.0 and k = 2.0. The 
size of each time-series was: T — 10^ with a time-step 
of M = 0.001 when r < 10 and T = lO" with a time- 
step was Ai = 0.01 when r > 10. When ,s = 1.0 and 
fc = 0.1 we used M = 1 in the region t < 1, A/ = 10 in 



4 



the region r G [1, 10] and M = 5 in the region r > 10. 
The size of each time-series was: T = 10^ with a time- 
step of At = 0.005 when t < 10 and T = 10" with 
a time-step was At — 0.01 when r > 10. The starting 
points of the simulated time-series were all the same with 
Xj{0) — 0.1 where j — 1, • • • , A/. The numerical simula- 
tions of the pdfs in the two bottom panels were obtained 
by using the same numerical data considered in the sim- 
ulation of the autocorrelation function when r > 10. In 
order to simulate the process in the y coordinate space, 
we started by simulating the process of Eq. ([T]) and com- 
puted y = fsix) for each simulated x value. 

By inspecting the two top panels of Fig. [2] it should 
be evident that the autocorrelation of the process defined 
by Eq. ([T]) with the coordinate transformation of Eq. Q 
can not be described by a single exponential, as in the 
case of the OU process. Indeed, by performing appropri- 
ate coordinate transformations one gets: 



^(A-^) / dyg{yT) sin{y){ll) 



=e-'=-/2Erf-i(l 



"'0 

-kz\ 



T 



4 



-1/2 



By using Lemma (6.1) in chapter 9 of Ref . (31) one can 
show that: 



T 



+00 2^ 

dyg{yT) sin(y) « — 

/o 

Therefore, for small energy values A we get: 



A^O 



(12) 



(13) 



This implies that for large time lags Rs{t) ~ 
exp(— ^t)t~^/^ as i — cx). We thus have a short-range 
correlated stochastic process whose autocorrelation func- 
tion is not exponential. 



B. A power-law correlated process 

Let us consider the stochastic process described by the 
following nonlinear Langevin equation [191 124j : 



X = h{x{t))+r{t) 

( -2yVbtan(v^x) if |x| < L, 
h{x) = <^ (14) 

[ (1 - ^/l + AVi)/x if \x\ > L. 

Vi^ L JVq tan(v/T^i) f 1 + i tan(v/VbL 



where L and Vq are real arbitrary constants and r(i) is 
a (5-correlated Gaussian noise term. It is straightforward 
to show that the stationary distribution of this process 
is: 

f iV7/cos(/K)x) if |a;| < L, 
Wi{x) = \ (15) 
[ Nil/x'^ if |a;| > L. 

a = y/l + AVi - 1 = 2L y/Vo tan(y/VoL) 



where Nj and Njj are real constants that can be obtained 
by imposing that Wi{x) is continuous and normalized to 
unity. By using the methodology of eigenfunction expan- 
sion [3] 2] it is possible to prove that the autocovariance 
function {x{t + r) x{t)) oc , where (3 — (a — 3)/2 thus 
showing that we are dealing with a power-law correlated 
stochastic process. In the range 3 < a < 5 the process is 
long-range correlated |19) . 

Let us now consider the coordinate transformation: 

2s Erf"' [1 - r{x)] \x\ > L 
2^ Erf"' [l~p{x)] \x\ < L 




(16) 



Ni 



(x-L) + 



2/V^L" cos(v%i)2 



[sm{2^/%x) ~ sm{2x/VoL)^ 



By using the Ito lemma, one can show that, starting from 
the process of Eq. ( 14 ), in the coordinate space y — fi (x) 
one gets a multiplicative stochastic process whose sta- 
tionary pdf is given by Eq. ([5]), with s as an additional 
arbitrary parameter. One can show that G(y) is con- 
tinuos along the real axis, although its first derivative is 
discontinuos in y — zLf{L). The drift coefficient H(y) is 
discontinuos iny = zLf{L). For large values of y the drift 
and diffusion coefficient behave as: 



1 



H{y) oc -sTrexp(- 
G{y) oc — ^exp(- 



y 



s{a — 1) 



2s(a- 1)' 



y 



According to the eigenfunction expansion methodology 
[31 Uj the autocorrelation function pi{t) of the process de- 
fined by Eq. ( 14 ) with the coordinate transformation of 
Eq. (16) is given by Eq. (l8| with fs{x) replaced by 



fi{x). Now the functions tpoixj and il^\{x) are the eigen- 
functions of the Schrodinger equation associated to the 
stochastic process of Eq. ( 14 ) and with potential: 



Vsix) 



-Vo if \x\ < L, 
Vi/x^ if \x\ > L, 



(17) 



In the top panels of Fig. [3] we show the results of nu- 
merical simulations of the autocorrelation function pi (t) 
performed for the case when s = 1.0, L = 1.0, and 
Vo = 0.98 (i.e. a = 3.05, top-left panel) and Vq = 1.16 
(i.e. a — 4.00, top- right panel). The solid (red) lines 
show the theoretical prediction obtained from Eq. ([s]) 
with fs{x) replaced by fi{x). The open circles show 
the results of the numerical simulations. By perform- 
ing a nonlinear fit (dashed blue line) , the autocorrelation 
function shows an asymptotic decay compatible with a 



5 




T > 10. When s = 1.0 and a = 4.00 we used M = 1 
in the region r < 1, Af = 20 in the region r > 1. The 
size of each time-series was: T = 10^ with a time-step of 
At = 0.005 when t < 10 and T = 10" with a time-step 
was At — 0.005 when r > 10. The starting points of the 
simulated time-series were all the same with Xj(0) = 0.1 
where j = 1, • • • , M. The numerical simulations of the 
pdfs in the two bottom panels were obtained by using the 
same numerical data considered in the simulation of the 
autocorrelation function when r > 10. In order to sim- 
ulate the process in the y coordinate space, we start by 
simulating the process of Eq. ( 14 ) and compute y = fi (x) 
for each simulated x value. 



The Doob Theorem 



FIG. 3: (Color online) The figure shows time-average numer- 
ical simulations performed according to Eq. ( |10[ ) for the case 
when when s — 1.0, and a — 3.05 (left), a — 4.00 (right). In 
the top panel we show the results for the autocorrelation func- 
tion. The solid (red) line shows the theoretical predictions 
obtained from Eq. ([8| with fsix) replaced by fi{x), while 
the open circles show the result of the numerical simulations. 
The dashed (blue) lines show a nonlinear fit performed on the 
large time region of the autocorrelation function. The bottom 
panel shows the stationary pdf of the process. Again the solid 
(red) line shows the theoretical prediction of Eq. ([sjl, while 
the open circles show the result of the numerical simulations. 
The parameters of the numerical simulations are given in the 
text. 



power-law r"''', with /?/ = 0.91 when a — 3.05 (top-left 
panel) and /?/ — 1.36 when a — 4.00 (top-right panel). 
The bottom panels of Fig. ([3 ) show the stationary pdf of 
the process. Again, the solid (red) line shows the theo- 
retical prediction of Eq. ([s]), while the open circles show 
the results of the numerical simulations for the case when 
a = 3.05 (bottom-left panel) and a = 4.00 (bottom-right 
panel) . 

By plotting the autocorrelation function pi (t) for L = 
1.0 and for different values of Vg, i.e. different values of 
a, and by performing a nonlinear fit, one finds that for 
large lag values pi{t) cx t^*^' where /3; shows a linear 
dependance from a described by = a/2 — 77 with 77 « 
0.61. By using the results of Ref. [23] one obtains that 
indeed the autocorrelation function decays like 
with logarithmic corrections. When a < 2(1 + 7]) the 
process is long-range correlated. 

Those shown in Fig. [3] are time-average numerical sim- 



10). The parame- 
ne autocorrelation 



ulations performed according to Eq 
ters of the numerical simulations of t 
function arc: AI = 1 in the region r < 1, Af = 10 in the 
region r G [1, 10] and Af = 20 in the region r > 10 when 
s = 1.0 and a = 3.05. The size of each time-series was: 
T = 10^ with a time-step of At = 0.005 when r < 10 
and T = 10^^ with a time-step was At — 0.005 when 



The existence of not exponentially correlated processes 
with Gaussian distribution does not contradict the Doob 
Theorem [3S]. In fact, such theorem deals with the case 
when the process admits a stationary pdf and a 2-point 
conditional transition probability that are both non sin- 
gular and Gaussian on the whole real axis. Indeed the 
Doob theorem applies to the Ornstein-Uhlenbeck (OU) 
process described by the following linear Langevin equa- 
tion [3]: 



i(t) = -jx{t) + DT{t) 
and whose conditional probability is given by [3]: 

P{x2,t2\xi,ti) 



(18) 



v/27ri:)2(l-e-27r) 



X exp( 



where r — t2 ~ ti. The two parameters D and 7 are 
related to the variance of the OU process by the relation 
s = DV7. 

A direct inspection of the conditional probabilities for 
the two processes of section II A and II B| shows that 
they are different from the one of Eq. (191. Indeed, a 



preliminary issue regards the way to compare the three 
processes considered above. We first notice that, by con- 
struction, the three processes have the same stationary 
pdfofEq. ^. This is a Gaussian distribution with vari- 
ance s. Hereafter we will always set s — I, which implies 
that D = in the OU process. As a result the OU pro- 
cess will be characterized by the single parameter 7, the 
process of section [II A| will be characterized by the single 
parameter k and the process of section [If B| will be char- 
acterized by the two parameters L and Vq or, by using Eq. 



( 15 1, the two parameters L and a. Without loss of gener- 



ality we can set L = 1. Secondly, we notice that the three 
parameters 7, k and a enter the autocorrelation function 
of the respective processes and therefore they determine 
how the processes loose or keep their memory. When 
comparing the three processes considered above one rea- 
sonable criterion for the choice of the parameters seems 
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the impose that they carry the same amount of memory, 
i.e. the integral of the autocorrelation function is the 
same. Of course this is one out of many possible choices, 
and we are also aware of the fact that it can not be ap- 
plied when considering long-range correlated processes. 
This choice however helps in presenting our results that 
are nevertheless based on analytical findings, as we will 
show in the next sections. 

In order to illustrate how the existence of not expo- 
nentially correlated processes with Gaussian distribution 
does not contradict the Doob Theorem, in Fig. |4] wc 
show a comparison between the conditional probability 
P{y,t\yo,0) of the stochastic process defined by Eq. ([T]) 
with the coordinate transformation of Eq. Q and the 
one predicted for the OU process. According to the above 
criterion, the parameters used for the process defined by 
Eq. ([!]) with the coordinate transformation of Eq. Q 
are s — 1.0, k = 2.0 and uq = 2.0. For the OU process 
we considered 7 = 2.039. The four panels show the be- 
haviour of the conditional probability P{y,t\yQ,0) with 
respect to y for four different values of time t: t = 0.1, 
t = 1.0, t = 5.0 and t = 10.0. The solid (red) lines are the 
theoretical predictions obtained from Eq. ( 19 1 and the 



open circles are the result of the numerical integrations 
performed according to 



P{y,t\yo,0)^Wg{y) + 

/fe2/4 



(20) 



dfs{x) %{xq) 

= fs\y) 



= /7'(yo) 



where "00(2^) and 'il)\{x) are the eigenfunctions of the 
Schrodinger equation with potential Vs{x)= — 
k5{x). Fig. [4] shows that for small values of time, the 
conditional probability for the process of section |II A| is 
quite different from the one of the OU process. As long as 
the time increases, the agreement gets better and better. 
In fact, it is expected that for large times both condi- 
tional probabilities must converge to the same pdf which 
is the one of Eq. ([5|. The bottom-right panel of Fig. |4] 
shows that at t = 10 the agreement is already good. 

In the next section we will analytically show that the 
way the conditional probability converges to the station- 
ary pdf is indeed different for the three processes. 



III. PERSISTENCE 

The existence of long-range interactions is often asso- 
ciated to the existence of persistencies in the time series 
of the stochastic process, i.e. to the fact that when a pro- 
cess reaches a position x then it spends some time around 
that position value. The conditional probability gives the 
possibility to investigate the occurrence of persistencies. 
In fact, we consider the conditional probability at equal 
spatial position P{x, t\x, 0) as a proxy of the persistence 
of the process. 




FIG. 4: (Color online) The figure shows a comparison between 
the conditional probability P{y,t\yo,Q) of the stochastic pro- 
cess defined by Eq. |[T]) with the coordinate transformation of 
Eq. Q and the one predicted for the OU process of Eq. ( |19[ l . 
Here, we consider the case when s = 1.0, k = 2.0 and j/o = 2.0. 
The four panels show the behaviour of the conditional proba- 
bility P{y, tjyo, 0) with respect to y for four different values of 
time t: t = 0.1, t = 1.0, t = 5.0 and t = 10.0. The open circles 
are the result of the numerical integrations of Eq. ( 20 1 and 



the solid (red) lines are the theoretical predictions obtained 
from Eq. ||19[), where we used 7 = 2.039. 



Starting from Eq 
P{y,t\y,0) for the processes 
tion [irei In Fig 



(201, we consider 



the quantity 
II A and sec- 



of section 
[5] we show the subtracted quantity 
Wg{y) for (i) the process of sec- 
II A (dashed line) with s — 1.0, k = 2.0, and {ii) 

1.0, 

a = 4.15 



V^jt) = Piy,t\y,0) 
tion 



section 

Vo = 1.183, i.e 



II B (dashed dotted line) with s — 1.0, L 

The solid line shows the same 



quantity for the OU process and it is obtained from Eq. 
(191 by simply imposing X2 = xi = t/o- The parameters 
used here are 7 = 2.039. In all the three cases we con- 
sidered y — yo — 2.0. The choice of these parameters 
ensures that the three processes have the same unitary 
variance and the same integral of the autocorrelation, i.e. 
they carry the same "amount of memory" . 

The figure shows how for small values of time the Vy (t) 
relative to the process of section [lIB| is higher than the 
other two. The functional behaviour seems simular to 
the behaviour which is typical of the OU case. In 

this range of time also the Vyit) relative to the process 
of section III Al seems to share the same behaviour as the 
OU process. The regime that however is interesting for 
our purposes is the one of large times. Here we observe 



that whilst the OU process and the process of section II A 
shows a very fast decay, the Vy (t) relative to the process 
of section |II B| shows a very slow decay that can be well 
described by a power-law with Py{t) « 

In Fig. |6] we show the function Vy (t) relative to the 
process of section II B for s = 1.0, L = 1.0 and different 
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on 

s 




2.95 
3.05 
3.50 
4.00 
4.50 
5.00 
5.50 

y=0.00271 l/x"0.8S2 
y=0.00214 l/x"0.912 
y=0.00U6 l/x"1.138 
y=0.00061 l/x"1.413 
y=0.00031 l/x«1.678 
y=0.000I5 l/x«1.937 
y=0.00006 l/x«2.191 



time t 



FIG. 5: (Color online) The figure shows a comparison between 
the subtracted conditional probability Vyit) = P{y,t\y,Q) — 
Wg[y) of (i)the OU process (black solid line) with parameters 
7 = 2.039, {a) the process of section |II A| (red dashed line) 
with parameters = 2.0, s = 1 and {iii) the process of section 

1.0, a = 

We show that whilst the OU process and 



II B (blue dash-dotted line) with parameters L 
4.00, s = 1.0 



the process of section II A shows a very fast decay, the Vy{t) 
relative to the process of section p^I B | shows a slow decay that 
can be well described by a power-law. 



FIG. 6: (Color online) The figure shows the subtracted condi- 
tional probability 'Py{t) = P{y, t\y, 0) — Wg{y) for the process 
of section p^I B | with parameters L — 1.0, s = 1.0 and different 
values of a, as shown in the legend. In the bottom part of the 
legend we report the result of a nonlinear fit performed in the 
range t £ [40,500]. These results show that Vyit) shows an 
asymptotic power behaviour Vy(t) ~ 1/i'', where /3 seem to 
depend linearly from a, according to the relation: /3 — a/2— ry. 
By using the results of asymptotic theory one can show that 
/? = (a-l)/2. 



values of Vq chosen in such a way that the parameter a 
assumes the values shown in the legend. In the bottom 
part of the legend we report the result of a nonlinear 
fit performed in the range t S [40,500]. These results 
show that Vy{t) shows an asymptotic power behaviour 
Vy{t) where /3 seem to depend linearly from a, 

according to the relation: /3 = a/2 — ry. 

Indeed, we are here interested in the large time be- 
haviour of the subtracted persistence function 'P{t) oc 
/q°° i^x{xf'^~^* , which is determined by the small A be- 
haviour of %j)\{xY [S]. One can show that in the limit 
when A — > 0: 



(21) 



Therefore the large time behaviour of V{t) is given by: 



Vy{t) « l/t 



(a-l)/2 



(22) 



which is in good agreement with the results of Fig. [6] 
Moreover, Eq. (22) confirms that the conditional proba- 



bility for such process is intrinsecally different from the 
one of Eq. ( 20 ) and therefore the Doob theorem can not 



be applied to this process. 

Analogously, o ne ca n also show that in the case of the 
process of section II A one gets ipxix) « (A— fc^/4)^/^ and 
therefore the subtracted persistence function Vy (t) shows 
a power-law decay with an exponential cut-off: Vy{t) w 

gain, this result analytically confirms the 
ones of section III CI 



It is worth mentioning that the integral dt Vy (t) al- 
ways exists for the OU process a nd th e process of section 
A For the process of section II B 'Py{t) is integrable 
only if a > 3, i.e. when the pocess is short-range cor- 
related. In this case, we can say that the conditional 
probability converges to the stationary distribution in a 
finite time. When a < 3, the process of section |II B| 
is long-range correlated, and also the above integral di- 
verges, thus suggesting that there is not a typical time in 
which the conditional probability converges to the sta- 
tionary distribution. Indeed, d,tVy{t) is related to 
the time spent by the process in a certain location. Such 
issue will be more quantitatively investigated in section 



Other processes with an infinite and unbounded 
set of timescales 



The power-law correlated stochastic processes men- 
tioned in the previous sections incorporates an infinite 
and unbounded set of timescales. This is due to the fact 
that the eigenvalues spectrum of the Schrodinger poten- 
tial of Eq. ( 17 1 is given by a continuum set of eigenvalues 



A > attached to the null eigenvalue associated to the 
ground state. As mentioned above, the power-law cor- 
related process shows a power-law decay in the persis- 
tence function Vy{t), as illustrated in Fig. [6j One might 
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therefore wander whether or not such power-law decay is 
merely due to the fact that the process incorporates an 
unbounded set of timescales. In general, any Schrodinger 
potential that asymptotically decays like would give 
a stochastic process with an infinite and unbounded set 
of time-scales. In the case /i = 1, we will show here that 
the corresponding persistence function does not necessar- 
ily display a power-law decay, . 

Let us now consider the additive stochastic process as- 
sociated to the Schrodinger potential [3] : 



Vs{x) 




(23) 



where L, Vo and Vi are real positive constants. This 
quantum potential is exactly solvable. In the region x > 
L the eigenf unctions are given by: 



^-0(2;) =Ao^ Ki{2 ^/Vi Vx) 



(24) 



Vl I r- 

i— -^,2,2iVXx) + 
2 



CaC/(1-«y -^,2,2i\/Ax))-hc. 



c. 



where Ki{-) is the K-Bessel function of order 1, U{-) is 
the confluent hypergeometric function and -F'i,i(-) is the 
Kummer confluent hypergeometric function. Usually the 
potentials arise in the context of quantum systems 
subject to a repulsive or attractive central Coulomb po- 
tential [26]. The Coulomb potential is then defined only 
for positive real values and the relevant eigenfunction 
must be regular at the origin. In our case, we are con- 
sidering a Coulomb potential extending over the whole 
real line although it is regularized in a; = by assuming 
the Vs{x) = — Vq when |a;| < L. This implies that the 
contribution associated to the confluent hypergeometric 
function U{-) must be taken into account. 

The normalization constant Aq is fixed by imposing 
that the ground state is normalized to unity. The other 
normalization constants A\ and Cx are fixed by imposing 
that the odd and even eigenf unctions fulfill the normal- 
ization condition J dx'ip\{x)'ip\>{x) = S{X — X'). The drift 
coefficient h{x) aymptotically decays like h{x) w l/y^- 

As mentioned above, the large time behaviour of the 
subtracted persistence function V(t) is determined by the 
small energy behaviour of '0a(2;)^ [31]. One can show that 
in the limit when A 0: 



and therefore the large time behaviour of Vxit) is: 



Vxit) 



^5/6 



(25) 



(26) 



The above result shows that the mere existence of an 
infinite and unbounded set of timescales is not enough 
to ensure an asymptotic power-law behaviour in Vx{t)- 



It is worth mentioning that such result is also valid for 
any coordinate transformation able to map the stochastic 



process associated to the quantum potential of Eq. ( 23 ) 
into a Gaussian distributed one. In passing, it is also 
possible to prove that the stochastic process associated 
to the quantum potential of Eq. ( 23 1 is not power-law 
correlated, even though it incorporates an infinite and 
unbounded set of timescales 1^. 



IV. MEAN FIRST PASSAGE TIME 



Naively speaking, the existence of power-law tails in 
Vy{t) indicates that the process with power-law autocor- 
relation function spends more time than the other two 
processes in a certain spatial position. This would imply 
that, starting from position x the power-law autocorre- 
lated process takes more time than the others to reach 
another position at distance A from x. In other words 
the process with power-law autocorrelation function is 
slower than the other two processes. Such an effect can 
be quantitatively measured by considering the mean First 
Passage Time [IZMH]- 

Specifically, we will consider the mean time (A) that 
is needed to reach for the first time position x±A starting 
from position x. This is obtained by solving the equation 



h{x) 



dx 



dx-^ 



(27) 



with boundary conditions (A) — when x = ± A. Here 
h{x) and g{x) are the drift and diffusion coefficient ap- 
pearing in the non linear Langevin equation associated 
to the considered stochastic process. 

In the case of the OU process of Eq. ( 18 ) the mean 
First Passage Time is given by: 



A^ 
2L>2 



H 



PFQ 



{l,l},{3/2,2},7 



A2 



2L>2 



(28) 



with D — ^7 s. 

For the stochastic Gaussian process of section II A the 



mean First Passage Time can be analytically computed 
by using the drift coefficient of Eq. ^ and the diffusion 
coefficient of Eq. ([7| into Eq. (27). For the case when 
X = one gets: 



To(A)--llog(l-Erf(A/\/2^; 



(29) 



k'^\ 1 - Erf (A/^ 



An alternative way to obtain such result is to analytically 
compute the mean First Passage Time for the stochastic 
process of Eq. ([T]) and then performing the substitution: 



A ^ log(l - Erf (A/y2s)) 



(30) 
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For large values of A the mean First Passage Time of Eq. 
(p9| gives: 



A 



(31) 



For the stochastic Gaussian process of section II B the 
mean First Passage Time can not be computed by ana- 
lytically solving Eq. (27). However, one can follow the 
alternative way mentioned above. Specifically, we first 
computed analytically the mean First Passage Time of 
the process of Eq. ([M] ) by using Eq. ( [27| . Subsequently, 
according to Eq. (|16p we performed the substitution: 



A^/7'(A) 



A^/7/(A) 



if A > L 
if A < L 



(32) 
(33) 



In the region A > fi{L) it is possible to obtain an ana- 
lytical expression showing that for large values of A one 
gets: 



ro(A)« (^e-')' 



A 



2s 



(34) 



Eq. (34) and Eq. (31) clearly show that for large values 
of A the process of section |IIB| has a mean First Passage 
Time larger than the process of section |II A[ This con- 
firms that the process that incorporates more time-scales 
is the slowest. 

For small values of A in all considered processes the 
mean First Passage Time grows quadratically. In fact, 
from Eq. (28) one easily gets that the OU process ad- 
mits To (A) « A^/(27s) while for the process of section 



II A from Eq. (29) one easily gets Tq{A) « k^/{-Ksk'^) in 



the limit when A — > 0. For the process of section |IIB| 
we do not have an analytical expression for To(A), since 
y = fii{x) is not analytically invertable although it is 
continuous and monotonic and therefore its inverse func- 
tion does exist. However, one can show that for small 
values of x one gets y = fii{x) ~ Kx. where: 



K = 



ALWo - a{l - 4L2Vb) - a cos(2LVK)) 



(35) 



By using the mean First Passage Time of the process of 
Eq. ( [T4| and by performing the substitution A i-^ A/K 
for small values of A one gets To(A) « A'^/{2K'^). 

In Fig. [Tjwe show the mean First Passage Time for the 
two processes of section II A (dashed line) , with param- 
eters s = 1.0, k = 2.00, and section II B| (dashed dotted 
line), with parameters s = 1.0, L = 1.0, Vq = 1.16, i.e. 
a = 4.00. For comparison, the solid line shows the same 
quantity for the OU process, with D = ^/j, 7 = 2.039. 
Again, the choice of these parameters ensures that the 
three processes have the same unitary variance and the 
same integral of the autocorrelation, i.e. they carry the 
same "amount of memory" . 

As expected, for large values of A the more time-scales 
the process includes, the larger the mean First Passage 




FIG. 7: (Color online) The figure shows a comparison be- 
tween the mean First Passage Time Tx{L) of (i)the OU pro- 
cess (black solid line) with parameters 7 = 2.039, {ii) the pro- 
cess of section |II A| (red dashed line) with parameters k = 2.0, 
s = 1 and (Hi) the process of section II B (blue dash-dotted 
line) with parameters L = 1.0, a = 4.00, s = 1.0. We show 
that for large values of A the more time-scales the process 
includes, the larger the mean First Passage Time. With this 
choice of the relevant parameters, although both processes of 
section [II A| and section [II B| are more persistent than the OU 
one, for small values of A they can be faster than OU. 



Time. Surprisingly, with this choice of the relevant pa- 
rameters, although both processes of section [Tl A[ and sec- 
tion [HB] are more persistent than the OU one, for small 
values of A they can be faster than OU. This suggests 
that when a process statistically spends more time in a 
given position, say a; = 0, the likelihood that it reaches 
nearby positions by chance is also enhanced. This is an- 
other way of looking at persistence. 



V. CONCLUSIONS 

The aim of the paper was to understand how the in- 
clusion of more and more time-scales into a stochastic 
process affects its conditional probability. We recall that, 
together with the stationary distribution, this is the other 
probability distribution that fully characterizes a station- 
ary Markovian process. 

In order to disentangle all possible effects due to the 
stationary distribution, we have considered processes 
with a Gaussian distribution and characterized by the 
presence of (z) one single time-scale, i.e. the OU process, 
(a) an infinite set of time-scales bounded from below that 
give rise to a short-range correlated process, see section 
II Al and 



{in) an infinite and unbounded set of time- 
scales that give rise to a power-law correlated process, 
see section Hi Bl 
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For these processes we have considered the equal po- 
sition conditional probability P{x,t\x^O) and the mean 
First Passage Time Tj;{A). The function P{x,t\x,0), or 
better Vxit) — P(x, t\x, 0) — Wg{x), tells us how the con- 
ditional probability converges to the stationary Gaussian 
distribution Wg{x). Moreover, P{x,t\x,0) can be con- 
sidered as a proxy of the persistence, i.e. the fact that 
when a process reaches a position x then it spends some 
time around that position value. The Mean First Passage 
Time tells us how much time needs a process to reach a 
position at distance A starting from position x. 

In section mil we have shown that the more time-scales 
the process includes, the larger the persistence. Specifi- 
cally, we have shown that the power-law correlated pro- 
cess of section|IIB shows a slow power-law decay of Vy (t) . 



correlated, the integral dt Vy (t) diverges. 



By contrast, the other two processes show a convergence 
dominated by an exponential cut-off. Moreover, in sec- 
tion |IIIA] we have shown that the existence of an infinite 
and unbounded set of timescales is a necessary but not 
sufficient condition for observing the power-law decay of 

Vyit). 

We have therefore given evidence that the existence of 
persistencies, as measured by Vyit), are neither necessar- 
ily related to existence of extreme events US] , nor to 
the existence of an infinite set of time-scales in the pro- 
cess. In the context of stationary Markovian processes 
the persistence seems to be associated to the power-law 
decay of the autocorrelation function, as for the process 
of section [ri B| This means that the weights by which each 
time-scale enters the process play a fundamental role. 
Moreover, it is worth mentioning that the power-law de- 
cay of Vy{t) can be observed for all values of a, i.e. when 
the process is either short-range and long-range power- 
law correlated. However, when the process is long-range 



In section IV we have shown that for large values of A 
the more time-scales the process includes, the larger the 
mean First Passage Time. Specifically, we have shown 



that the process of section pTB shows a power-law growth 
while the process of secton II A shows a linear growth in 
terms of the (ze^ ) function, with z = A/\/2s. This 
indicates that, starting from position x = 0, locations 
that are very distant are reached with increasing diffi- 
culty. That raises the problem of the effective capacity 
of a power-law correlated process to span all the real 
axis in a finite time. Surprisingly, for small values of A, 
at least for some choices of the relevant parameters, pro- 
cesses that include more time-scales show a smaller mean 
First Passage Time. This suggests that when a process 
statistically spends more time in a given position, say 
X — 0, the likelihood that it reaches nearby positions by 
chance is also enhanced. 

Further work is needed to understand how the above 
results are related to the ergodicity properties [34 of the 
process of section |IIB[ Specifically, a issue worth of fur- 
ther investigation is how (z) the divergence of dt Vy (t) 
when a < 3 and (ii) the power-law growth of the mean 
First Passage Time affect the fact that the process is er- 
godic or not. 
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